F/6  */l 


AD-A113  1*2  AIR  FORCE  GEOPHYSICS  LAB  HANSCOM  AFB  MA 

A  6ROUNO  TRUTH  ANALYSIS  OF  OMSP  WATER  VAPOR  RADIANCES. (U> 
NOV  SI  F  R  VALOVCIN 
UNCLASSIFIED  AFGL-TR-81-0323 


MICROCOPY  RESOLUTION  TEST  CHART 


Unclassified 


?0»  (Coot) 

~4;  1  ratio,  is  found  in  all  the  DMSP  SSH  H2O  channels.  These  results  are  in 
agreement  with  MeClatchey's  1876  results  in  his  analysis  of  the  DMSP  15  am 
CO2  sounder  channels.  ^ 

The  data  sets  comparison  were  divided  into  three  latitude  belts,  that  is. 
Tropical,  Mid-latitude  and  Arctic.  In  turn,  the  discrepancies  between  cal¬ 
culated  and  measured  radiances  appear  to  be  latitudinally  dependent.  Smaller 
discrepancies  are  found  in  the  Tropica  and  the  larger  discrepancies  are  found 
in  the  Arctic  latitude  belt.  Also,  it  appears  that  the  DMSP  SSH  HoO  channels 
cannot  discriminate  between  low  cloud  contamination  and  clear  column  condi¬ 
tions. 
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A  Ground  Truth  Analysis  of  DMSP 
Water  Vapor  Radiances 

1.  INTRODUCTION 

The  Defense  Meteorological  Satellite  Program  (DMSP)  Block  5D  Satellite  carries 
a  Special  Meteorological  Sensor  H  (SSH)  package.  An  excellent  review  of  the  optical 
subsystems  and  the  spectral  characteristics  of  this  SSH  package  is  provided  by 
Nichols.  Infrared  energy  is  measured  in  16  spectral  bands  by  the  SSH  package: 
six  spectral  channels  are  located  in  the  15  Jim  carbon  dioxide  bands,  eight  in  the 
18  to  30  pm  rotational  water  vapor  band,  one  in  the  9.  6  pm  ozone  band,  and  one 
in  the  atmospheric  window  near  12  pm.  In  this  study,  the  emphasis  is  placed  on 
the  eight  DMSP  water  vapor  SSH  channels.  The  various  channel  characteristics 
are  listed  in  Table  1. 

Originally,  the  purpose  of  this  study  was:  (1)  to  evaluate  the  DMSP  multi¬ 
channel  water  vapor  radiances  for  informational  content,  and  (2)  to  recommend 
some  operational  techniques  whereby  moisture  parameters  may  be  derived  directly 
from  DMSP  water  vapor  radiance  measurements.  A  preliminary  evaluation  of  the 
DMSP  water  vapor  radiances  was  made  in  a  previous  report2  concerning  the 

(Received  for  publication  2  November  1981) 

1.  Nichols,  D.  A.  (1975)  DMSP  Block  5D  special  meteorological  sensor  H,  optical 

subsystem.  Opt.  Eng.  14:284-288. 

2.  Valovcin,  F.  R.  (1980)  DMSP  Water  Vapor  Radiances— A  Preliminary 

Evaluation,  AFGL-TR-80-0312,  AD  AOSSSOT. 


informational  content.  Before  a  reliable  operational  technique  on  deriving  moisture 
parameters  could  be  recommended,  the  Forward  Problem  had  to  be  investigated. 

By  definition,  the  Forward  Problem  is  the  matching  of  satellite  observed  or 
measured  and  the  classically-calculated  radiances  based  on  coincident  in  time  and 
space  radiosonde  data.  McClatchey  in  studying  the  Forward  Problem  reported  a 
systematic  discrepancy,  that  is,  calculations,  in  general,  exceed  the  measured  or 
observed  radiances  in  the  DMSP  IS  pm  carbon  dioxide  bands.  As  this  study  pro¬ 
gressed,  it  became  quite  apparent  that  the  purpose  would  change  from  recommending 
some  operational  techniques  on  deriving  moisture  parameters  to  conducting  a 
Ground  Truth  Analysis  of  the  DMSP  water  vapor  radiances  similar  to  McClatchey' s 
study  of  the  DMSP  15  pm  carbon  dioxide  band.  Thus  in  this  report,  a  comparison 
program  was  carried  out  for  approximately  70  sets  of  DMSP  SSH  water  vapor 
measurements  with  calculated  water  vapor  radiances.  Both  clear  and  cloud  con¬ 
taminated  columns  as  a  function  of  three  latitude  belts,  that  is.  Tropical,  Mid¬ 
latitude  and  Arctic  were  used  in  this  study.  The  DMSP  water  vapor  radiance  sets 
analyzed  were  from  Flight  II  SSH  aboard  spacecraft  WX  13536  and  on  a  few  occa¬ 
sions  from  Flight  III  SSH  aboard  spacecraft  WX  14537.  Flight  II  SSH  and 
Flight  HI  SSH  were  launched  in  July  1977  and  April  1978  respectively. 


Table  1.  DMSP  Water  Vapor  SSH  Channel 
Characteristics* 


Band 

pm 

Center  _  j 

cm 

Width 

cm-1 

NESR** 

FI 

28.  2 

355.0 

15.0 

0.  25 

F2 

25.  2 

397.5 

10.0 

0.  16 

F3 

23.8 

420.  0 

20.0 

0.  12 

F4 

22.7 

441.5 

18.0 

0.09 

F5 

18.7 

535.0 

16.  0 

0.  15 

F6 

24.5 

408.5 

12.0 

0.  14 

F7 

26.7 

374.0 

12.0 

0.  18 

F8 

28.3 

353.5 

12.  0 

0.33 

After  Nichols1 

>jc* 

NESR:  Noise  Equivalent  Spectral  Radiance  in 
mW/mz  sr  cm'* 


3.  McClatchey,  R.A.  (1976)  Satellite  Temperature  Sounding  of  the  Atmosphi 
Ground  Truth  Analysis,  APGL-TR-7o-fl279,  Ab  A038236. 


ere: 


10 


H 


^ija'ii  mi 


2.  DMSP  WATER  VAPOR  RADIANCE  DATA 

The  Air  Force  Global  Weather  Central  (AFGWC)  forwarded  to  the  Air  Force 
Geophysics  Laboratory  (AFGL)  printouts  of  world-wide  SSH  spectral  radiance  data 
for  the  period  00Z  30  January  1979  to  12Z  18  October  1979.  The  printouts  con¬ 
tained  information  on  date -time,  latitude  and  longitude  of  the  scan  spot.  Zenith 
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angle,  3D  NEPH  and  data  base  parameters  along  with  the  spectral  radiance 
measurements  for  the  16  DMSP  spectral  channels.  The  scan  spot  is  the  first  cloud- 
free  radiance  set  available  within  a  3°  X  3°  area  as  selected  by  AFGWC's  HPKG 
software  (based  on  the  3D  NEPH  files).  Only  the  eight  water  vapor  channels  and 
the  one  located  in  the  atmospheric  window  near  12  pm  were  used  in  this  study. 

The  DMSP  SSH  sounder  scan  pattern  is  accomplished  in  25  4 -degree  incremental 
steps  from  -48°  to  +48°  of  nadir  across  the  orbit  track.  The  range  of  zenith  angles 
for  the  25  scan  spots  are  +57®  to  -57°.  Zenith  angles  are  designated  positive  on 
the  left  and  negative  on  the  right  of  the  scan  swath.  At  nadir  the  scan  spot  on  the 
earth's  surface  is  approximately  21  nmi  (39  km)  and  increase  to  approximately 
35  by  25  nmi  (65  X  46  km)  at  a  zenith  angle  of  +/-32°.  The  analysis  was  limited 
to  zenith  angles  of  +/-32°  or  approximately  +/- 245  nmi  (454  km)  from  nadir. 


3.  GROUND  TRUTH  DATA 

The  USAF  Environmental  Technical  Application  Center  (ETAC)  provided  the 
ground  truth  data  for  the  various  stations  world -wide.  A  "ground  truth"  site  is  de¬ 
fined  as  a  radiosonde  observation  being  within  +/-3  hrs  and  100  nmi  (185  km)  of  a 
DMSP  satellite  scan  spot.  Atmospheric  conditions  such  as  clear  or  cloud  con¬ 
taminated  were  determined  through  surface  observations,  AFGWC  three-dimensional 
nephanalysis  (3D  NEPH)  and  the  radiosonde  observations.  Imagery  was  not  gener¬ 
ally  used  because  of  the  period  involved  (January -October  1979)  and  the  final  selec¬ 
tion  of  the  "ground  truth"  site  world-wide.  The  radiosonde  data  were  received  in 
a  LOWTRAN  format  with  the  following  parameters:  Height  (km).  Pressure  (mb). 
Temperature  (®C),  Dewpoint  Temperature  (®C),  Relative  Humidity  (%),  and 
Absolute  Humidity  (g/m  ).  The  values  of  absolute  or  relative  humidity  above 
300  mb  are  climatologically  modeled  by  ETAC  since  the  radiosonde  usually  does 
not  report  humidity  above  this  level.  An  example  of  the  detailed  atmospheric  pro¬ 
file  of  temperature,  water  vapor,  and  ozone  as  a  function  of  height  and  pressure 
that  is  used  in  the  transmittance  calculations  is  shown  in  Table  2.  The  concentration 

4.  Coburn,  A.  R.  (1970)  Three  Dimensional  Analysis,  AFGWC,  Offut  AFB,  NE, 

AFGWC  TM-70-9. 
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of  water  vapor  was  calculated  as  described  by  Selby  and  McClatchey.  The  ozone 

data  were  introduced  from  climatological  models  of  ozone  distributions.  The 

2 

units  of  water  vapor  and  ozone  are  in  molecules/cm  and  are  integrated  values 
from  the  top  of  the  atmosphere  to  the  surface. 


Table  2.  Atmospheric  Profile  and  Composition  Derived  From 
Radiosonde  Data— West  Palm  Beach,  Florida,  72203-790403 


HT 

(km) 

Pressure 

(mb) 

Temp 

<K) 

Water 

Vapor 

2 

(Molecules /cm  ) 

Ozone 

2 

(Molecules /cm  ) 

31.  3 

10.  0 

237.  9 

1.  067E+19 

3.  451E+16 

29.  2 

13.  5 

232.  5 

1.  155E+20 

6.  160E+17 

28.  7 

14.  5 

235.  3 

1.  372E+20 

8.  015E+17 

26.6 

20.  0 

228.  1 

2.  106E+20 

1.  815E+18 

25.6 

'  23. 0 

223.  1 

2.  230E+20 

2.  259E+18 

23.  9 

30.0 

221.  7 

2.  377E+20 

3.  117E+18 

20.  7 

50.0 

211.  9 

2.  496E+20 

4.  764E+18 

19.4 

62.0 

200.  3 

2.  517E+20 

5.  439E+18 

18.  7 

70.0 

200.  9 

2.  529E+20 

5.  819E+18 

17.8 

81.0 

200.  3 

2.  544E+20 

6.  287E+18 

17.7 

82.0 

198.  1 

2.  545E+20 

6.  327E+18 

16.  9 

95.0 

203.  9 

2.  561E+20 

6.  78813+18 

16.6 

100.  0 

201.  9 

2.  569E+20 

6.  947E+18 

16.4 

102.0 

201.  3 

2.  571E+20 

7.008E+18 

14.6 

138.0 

202.  3 

2.  661E+20 

7. 859E+18 

14.  1 

150.0 

203.  9 

2.  625E+20 

8.  099E+18 

13.6 

165.  0 

205.  9 

2.  646E+20 

8. 329E+18 

12.4 

200.0 

215.7 

2.  764E+20 

8.  95712+18 

10.9 

250.0 

227.  9 

3.  396E+20 

9.  639E+18 

10.2 

278.0 

236.  3 

4.  905E+20 

9.  968E+18 

9.7 

300.  0 

237.  9 

6.  983E+20 

1.  020E+19 

7.6 

400.0 

254.7 

1.  266E+21 

1. 110E+19 

5.  9 

500.  0 

266.3 

1.  949E+21 

1.  282E+19 

4.7 

581.0 

275.  6 

3.  294E+21 

1.  231E+19 

3.  2 

760.0 

284.  0 

6.  999E+21 

1. 291E+19 

2.  2 

791.  0 

289.  6 

1.  125E+22 

1. 330E+19 

1.  9 

811.0 

289.  6 

1.  233E+22 

1.  338E+19 

1.9 

819.0 

283.0 

1.  336E+22 

1.341E+19 

1.5 

850.0 

285.  8 

2.  164E+22 

1.  352E+19 

1.  3 

876.0 

287.0 

2.  877E+22 

1.  360E+19 

0.7 

944.0 

290.4 

5. 222E+22 

1.  379E+19 

0.2 

1000.  0 

295.  2 

7. 703E+22 

1.  393E+19 

0.0 

1018.  0 

295.4 

8.475E+22 

1.  396E+19 

These  are  integrated  values 
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4.  ATMOSPHERIC  WATER  VAPOR  TRANSMITTANCE 


Atmospheric  transmittance  calculations  for  the  DMSP  SSH  water  vapor  channels 
for  the  approximately  70  selected  cases  were  computed  using  the  McClatchey  trans¬ 
mittance  program  developed  at  AFGL.  In  this  program,  line  by  line  calculations 

3 

are  performed  using  the  following  equations  from  McClatchey. 

1.0 

I.,  =  [  /  /  B<„.V)dr  d,+  /  f(t')B(j/,  T„)  dtz)  /  /  Uu)6u  (1) 

Ag  Tg  Ay  a  /  Al» 

which  is  the  solution  of  the  radiative  transfer  equation,  where 

2  -1 

I  is  the  radiant  intensity  in  mW/m  sr  cm  , 

B(p, T)  is  the  Planck  blackbody  function, 

T  is  the  atmospheric  temperature  and  Tg  is  the  surface  temperature, 

T  is  the  transmittance  of  the  atmosphere  from  the  altitude  associated 

with  the  pressure  level,  p,  to  the  top  of  the  atmosphere, 
v  is  the  frequency  (given  here  in  cm  1 ),  and 

f(y)  is  the  DMSP  SSH  instrument  filter  function. 


McClatchey  assumes  that  B(y.T)  is  relatively  constant  over  the  width  of  a  filter 
function  (10  or  20  cm  1  wide),  and  he  writes  ft]  (p)  as  an  independent  variable  in¬ 
stead  of  r,  and  obtains  Eq.  (2),  where  the  quantity,  dr/dfT)  p.  now  becomes  a 
weighting  function  that  can  be  interpreted  as  defining  the  atmospheric  layer  pri¬ 
marily  responsible  for  the  upwelling  emission  in  the  spectral  interval.  Ay, 


^  *  /  B(];'T)awpF  d£l>p  +  B(T;-  V 

pg 


-  _  f(v)r  (v) dv) 

A  Uv)du  ‘ 


Monochromatic  transmittances  were  computed  over  the  various  eight  DMSP  water 
vapor  filter  functions  for  the  appropriate  radiosonde  observation  taking  tempera¬ 
ture,  water  vapor  and  ozone  distributions  as  a  function  of  pressure  and  height  as 
shown  in  Table  2.  These  monochromatic  transmittances  were  then  weighted  by  the 
appropriate  filter  function  as  indicated  in  Eq.  (3)  in  order  to  generate  the  appro¬ 
priate  averaged  transmittance.  The  AFCRL  Atmospheric  Absorption  Line 

O 

Parameters  Compilation  was  used  for  all  water  vapor  absorption  lines  between 

6.  McClatchey,  R.A.,  Benedict,  W.S. ,  Clough,  S.  A.,  Burch,  D.  E. ,  Calfee,  R.  P. , 
Fox,  K. ,  Rothman,  L.S.,  and  Garing,  J.  S.  (1973)  AFCRL  Atmospheric 
Absorption  Line  Parameters  Compilation,  AFCRL-tR-V3-c)096,  Ad  3732904. 
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320-600  cm  .  The  calculation  was  based  on  the  April  1979  data  tape.  The 

Lorentz  line  shape  was  used  throughout  with  a  line-wing  modification  of  the  water 

7 

continuum  proposed  by  Burch  and  Gryvnak.  The  simple  Lorentz  line  shape  is  used 
for  1 1/- i/o |  ^  20  cm"1  and  varies  linearly  from  unity  at  1 -  t-o  |  =  20  cm-1  to  zero 
at  1 1/  -  k> I  =  30  cm  .  Also  from  this  same  report  by  Burch  and  Gryvnak,  the 
appropriate  values  of  the  temperature  coefficient  for  the  water  vapor  continuum 
were  used  and  are  listed  in  Table  3. 


Table  3.  Temperature  Coefficient  for 
Empirical  Water  Vapor  Continuum* 


V  . 

(cm"*) 

b 

(K) 

V 

(cm-1) 

b 

(K) 

320. 

1197 

330. 

1206 

1214 

350. 

1223 

600. 

1442 

360. 

1231 

610. 

1454 

370. 

1240 

620. 

1467 

380. 

1248 

630. 

1481 

390. 

1257 

640. 

1495 

400. 

1265 

650. 

1510 

410. 

1273 

660. 

1525 

420. 

1282 

670. 

1542 

430. 

1290 

680. 

1560 

440. 

1299 

690. 

1579 

450. 

1307 

700. 

1597 

460. 

1316 

710. 

1615 

47  0. 

1324 

7  20. 

1633 

480. 

1333 

730. 

1652 

49C. 

1341 

740. 

1670 

500. 

1349 

750. 

1688 

510. 

1358 

760. 

1706 

520. 

1366 

770. 

1725 

530. 

1375 

780. 

1743 

540. 

1383 

790. 

1761 

550. 

1392 

800. 

1779 

560. 

1400 

810. 

1798 

570. 

1409 

820. 

1816 

580. 

1419 

590. 

1430 

ac,  (0)  =  ac„  (296  K)  exp  (*■  - 
s  s 

$ 

After  Burch  and  Gryvnak 

Burch,  D.  E.,  and  Gryvnak,  D.  A.  (1979)  Method  of  Calculating  H?0  Trans 
mittance  Between  333  and  633  cm'*),  winai  nsssSTTTTITT  — 


S.  DMSP  WATER  VAPOR  FILTER  FUNCTIONS 


In  Appendix  filter  transmission  curves  and  the  digitized  filter  functions  for 
the  eight' DMSP  water  vapor  channels  are  provided.  These  curves  and  filter  func¬ 
tions  are  valid  for  the  Flight  II  SSH  package  aboard  spacecraft  WX  13536  launched 
in  July  1977.  The  listed  frequencies  for  each  channel  are  normally  the  central 
frequency  for  each  filter.  The  frequency  steps  for  each  channel  were  0.5  cm  *. 
These  curves  were  used  to  calculate  the  water  vapor  transmittances  and  the  result¬ 
ing  weighting  functions  according  to  Eqs.  (2)  and  (3). 

6.  COMPARISONS  BETWEEN  DMSP  WATER  VAPOR 

MEASUREMENTS  AND  CALCULATED  RADIANCES 

The  DMSP  water  vapor  measurements  were  provided  by  AFGWC  from  their 

HPKG*RADISAVE  HPKG  files  for  the  period  February  to  October  1979.  The 

radiances  set  of  DMSP  SSH  data  is  the  first  cloud-free  set  available  within  a 

3°  X  3°  area  as  selected  by  AFGWC's  HPKG  software  (based  on  3D  NEPH  files). 

o 

As  described  in  an  earlier  report,  this  file  also  contained  cloud  contaminated  scan 
spots.  The  radiosonde  observation  which  was  used  as  the  "ground  truth"  site  was 
obtained  from  ETAC  in  a  LOWTRAN  format.  The  transmittance  calculations  were 
based  on  the  temperature  and  humidity  profiles  as  a  function  of  pressure  from 
these  radiosonde  observations.  The  major  criteria  used  in  the  colocated  DMSP 
water  vapor  data  and  the  "ground  truth"  site  were  as  follows:  (1)  radiosonde  sta¬ 
tion  located  within  100  nmi  (185  km)  of  a  DMSP  scan  spot,  (2)  radiosonde  observa¬ 
tions  taken  within  +/ -3  hr,  and  (3)  zenith  angle  of  the  scan  spot  restricted  to 
+/- 32°. 

The  comparison  between  the  DMSP  water  vapor  measurements  and  the  calcu¬ 
lated  radiances  were  divided  into  three  latitude  belts;  that  is,  Tropical  (26°S-26°N), 
Mid-latitude  (26',-62°N)  and  Arctic  (62°-90°N).  In  addition,  the  ground  truth 
analysis  considered  both  the  clear  and  cloud  contaminated  scan  spots.  For  illus¬ 
tration  purposes.  Figures  1,  2,  and  3  depict  the  theoretical  weighting  function 
curves  for  the  eight  DMSP  water  vapor  channels  for  a  Tropical,  Standard  and 
Sub-arctic  atmosphere  respectively.  The  pressure,  temperature  and  water  vapor 

g 

concentration  for  these  three  models  are  thosed  used  by  McClatchey  at  al.  On 
the  left  side  of  these  figures,  the  most  opaque  water  vapor  channels  (F8-353, 


8.  McClatchey,  R.A.,  Fenn,  R.W. ,  Selby,  J.  E.  A.,  Voltz,  F.  E. ,  and 

Garing,  J.S.  ( 1972)  Optical  Properties  of  the  Atmosphere  (Third  Edition), 
AFCRL-72-0497,  AD  A753075. 


Fl-355,  F7-374  and  F2-397  cm"*)  are  shown.  The  Less  opaque  or  more  trans¬ 
parent  water  vapor  channels  (F3-420,  F4-441,  F6-408,  and  F5-535  cm"*)  are 
shown  on  the  right.  The  most  opaque  DMSP  water  vapor  channel  is  F8-353  cm  * 
and  the  most  transparent  is  F5-535  cm'1.  The  maximum  value  of  the  weighting 
function  for  each  channel  represents  the  approximate  location  in  the  atmosphere 
from  which  the  major  portion  of  the  energy  is  received  at  the  satellite  sensor. 

Table  4  lists  the  location  in  mb  of  the  maximum  value  of  the  weighting  function 
for  the  three  model  atmospheres. 

6.1  Gear  Column  Water  Vapor  Radiance  Compariaon 

The  colocated  cases  selected  for  comparison  purposes  designated  clear  are 
listed  in  Tables  5,  6,  and  7  for  Tropical,  Mid-latitude,  and  Arctic  latitude  belts 
respectively.  The  40  clear  cases -Tropical  (10),  Mid-latitude  (21),  and  Arctic  (9) 
were  selected  based  on  surfac  .■  reports,  3D  NEPH  and  the  radiosonde  humidity 
profile.  The  tables  list  the  stations,  location,  date,  zenith  angle,  distance  be¬ 
tween  the  satellite  scan  spot  and  the  radiosonde  or  "ground  truth"  station  along 
with  the  satellite  measured  and  the  calculated  water  vapor  radiances  in 
mW/m  sr  cm  for  the  eight  DMSP  water  vapor  channels.  The  calculated  radiances 
were  computed  from  Eqs.  (2)  and  (3)  using  the  radiosonde  of  temperature  and 
humidity  observations  coincident  in  space  and  time  as  "ground  truth". 

Figures  4,  5,  and  6  show  the  results  of  the  comparison  in  graphical  form.  As 
can  be  seen  from  both  the  tables  and  the  figures,  the  calculated  water  vapor  radi¬ 
ances  are  greater  than  those  measured  by  the  satellite  in  the  majority  of  cases. 

The  percentage  deviations  in  the  mean  are  listed  in  the  tables  and  are  defined  as 
(measured  -  calculated)  ~  measured  radiance  in  percent.  A  positive  (negative) 
deviation  indicates  that  the  measured  is  greater  (less)  than  the  calculated  radiance. 

In  the  mean,  the  smallest  percentage  deviations  are  found  in  the  Tropics  and  the 
largest  in  the  Arctic  latitude  belt.  The  range  of  mean  percentage  deviations  for  the 
eight  DMSP  channels  are  +0.  2  percent  to  -4.  5  percent  for  the  three  latitude  belts 
in  the  clear  column  comparisons.  The  ratio  of  negative  to  positive  percentage  de¬ 
viations  for  the  various  latitude  belts  are:  Tropic  52:28,  Mid-latitude  132:36,  and 
Arctic  70:2  for  a  total  of  253  negative  vs  67  positive  percentage  deviations.  Thus, 
the  calculated  as  indicated  by  the  negative  deviation  exceed  the  satellite  measured 
water  vapor  radiances  by  almost  a  4:1  ratio.  In  addition,  the  calculated  exceed  the 
measured  water  vapor  radiances  in  excess  of  +/-2.0  mW/m  sr  cm  by  almost  a 
3:1  ratio.  The  correlation  between  measured  and  calculated  water  vapor  radiances 
are  higher  for  the  most  transparent  channels,  that  is,  F5-535  and  F6-408  cm'1  and 
lower  for  the  most  opaque  water  vapor  channels,  that  is,  F8-353  and  Fl-355  cm'1. 
The  lower  correlations  for  the  most  opaque  channels  may  be  due  to  the  use  of  the 
climatological  humidity  values  above  300  mb.  Also,  it  should  be  noted  the  maximum 
percentage  deviations  for  all  three  latitude  belts  are  found  in  the  DMSP  water  vapor 
channel  F3-420  cm'1. 
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Table  4.  Location  of  the  Maximum  Value  of  the  Weighting 
Function  for  the  DMSP  Water  Vapor  Channels 


DMSP  Channel 

Tropical 

(mb) 

Standard 

(mb) 

Arctic 

(mb) 

FI -355  cm'1 

355 

385 

480 

F2-397 

405 

580 

555 

F3-420 

525 

580 

635 

F4-441 

525 

660 

730 

F5-535 

760 

845 

835 

F6-408 

595 

660 

835 

F7-374 

460 

505 

635 

F8-535 

355 

385 

480 

Table  5.  Comparison  of  Measured  With  Calculated  Radiances— Tropical 
(26S-26N)— Clear 


Location 

Date 

9 

Distance 

FI 

355 

F2 

397 

F3 

420 

F4 

441 

F5 

535 

F6 

408 

F7 

374 

F8 

353 

Measured 

Calculated 

Jiddah,  SA 

21.  5N  39.  2E 

790406 

22.74 

99  nmi 

79.  0 
81.5 

92.  9 
94.5 

100.0 

104.9 

107.5 

110.6 

130. 1 

129.2 

105.6 

106.0 

86.6 

90.  5 

74.8 

77.9 

Measured 

Calculated 

Jiddah,  SA 
'2I.5W  39.  2E 

790406 

18.  16 

66  nmi 

79.0 

81.7 

92.4 

94.7 

100.0 
105.  1 

107.6 

110.8 

129.4 

129.4 

105.6 

108.2 

86.5 

90.6 

74.5 

78.1 

Measured 

Calculated 

Chiang  Mai.  TH 
18.78N  98.98E 

79D409 

18.  16 

90  nmi 

75.  1 
75.7 

88.3 

87.4 

94.6 

96.6 

101.8 

102.0 

120.  1 
119.7 

100.6 

100.5 

82.3 

83.  1 

71.4 

72.7 

Measured 

Calculated 

Ascension  Island 
7.97S  14. 4W 

7 90329 

9.  05 

99  nmi 

82.  7 
77.4 

94.3 

89.8 

102.0 

99.6 

108.4 

105.0 

124.4 
123.  1 

105.5 

103.5 

89.  1 

85.  4 

78.6 

74.2 

Measured 

Calculated 

Ascension  Island 

7.  97S  14. 4W 

790507 

22.8 

77  nmi 

82.  1 
80.9 

94.6 

92.  8 

100.9 

101.7 

107.7 

106.6 

122.1 

121.9 

105.0 

104.3 

89.  1 
88.8 

82.8 

77.8 

Measured 

Calculated 

Coolidge  Field.  AT 
17.  12N  61.78W 

790202 

-9.  05 

52  nmi 

84.7 

81.  2 

95.4 

93.0 

103.8 
101. 1 

110.9 

107.0 

126.7 
126. 1 

107.6 

105.8 

92.2 

88.  1 

79.6 

77.2 

Measured 

Calculated 

TRUK  KA 

7.47N  151. 8E 

790603 

-18.04 

31  nmi 

78.0 

76.7 

88.5 

89.  1 

94.5 

98.6 

101.3 
104.  1 

119.  1 
125.0 

99,3 

103.2 

82.  9 

84.  1 

73.5 

73.7 

Measured 

Calculated 

MT  ISA  MO  AU 
20.67S  139. 48E 

790328 

-22.  86 

92  nmi 

75.0 

78.5 

89.5 

90.5 

96.2 
100.  1 

104.2 

105.3 

127.4 

124.2 

103.0 

103.4 

82.4 

86.5 

70.3 

75.4 

Measured 

Calculated 

MT  ISA  MO  AU 
20.67S  139. 48E 

790330 

32.  2 

32  nmi 

76.  3 
76.9 

88.5 

88.6 

94.6 

97.7 

101.2 

102.8 

119.2 
120. 1 

99.5 

101.2 

82.8 

84.4 

72.0 

74.0 

Measured 

Calculated 

Giles  MO  AU 

25.03S  129.  90S 

790401 

4.52 

91  nmi 

73.0 

75.  8 

84.3 

86.8 

90.0 

95.4 

96.3 

100.3 

113.6 
119.  1 

94.8 

98.9 

78.7 

62.8 

69.2 

73.2 

Measured  Mean 
STD  DEV 

78.3 

3.9 

90.  9 

3.6 

97.7 

4.3 

104.7 

4.5 

123.2 

5.3 

102.7 

4.0 

85.3 

4. 1 

74,7 

4.4 

Calculated  Mean 
STD  DEV 

78.6 

2.5 

90.7 

2.9 

100. 1 
3.2 

105.5 

3.4 

123.8 

3.7 

103.7 

3.0 

86.4 

2.9 

75.4 

2. 1 

(Meas-Calc)  „ 

-0.4 

•fO.  2 

-1.4 

Meas 

Correlation 

0.  662 

0.780 

0.724 

0.784 

0.850 

0.833 

0. 657 

-0.  594 
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Table  6.  Comparison  of  Measured  With  Calculated  Radiances— Mid-latitude 
(26-62N)  and  (26-62S)-Clear 


FI 

F2 

F3 

F4 

F5 

F6 

F7 

F8 

■ 

Location 

Date 

6 

Diatance  355 

397 

420 

441 

535 

408 

374 

353 

Lajes,  Azores 
38.73N  27.08W 

790711 

8.  99* 

63  nmi 

78.5 

83.4 

88.4 

95.8 

95.6 

105.2 

102.5 

110.4 

120.0 

125.7 

99.7 

107.8 

64.8 

91.6 

73.2 

80.3 

Grlefswald.  Denmark 

791001 

9.05* 

18  nmi 

72.  2 

81.7 

89.  1 

95.6 

111.9 

93.6 

78.  7 

67.6 

54.  ION  13.  38E 

71.7 

82.4 

90.3 

95.0 

108.8 

94.2 

78.  2 

69.2 

OMSK,  USSR 

54.  93N  73.40E 

790314 

0* 

59  nmi 

68.3 

70.7 

77.9 

60.0 

82.3 

86.5 

87.  1 
90.3 

92.4 

97.8 

86.5 

89.1 

73.2 

7  8.8 

65.8 

68.5 

OMSK.  USSR 

54.  93N  73.40E 

780505 

27. 27* 

45  nmi 

75.  1 
72.7 

65.4 

82.  5 

90.6 

89.7 

95.9 

93.9 

108.2 

105.5 

94.2 

92.5 

80.6 

79.0 

71.8 

70.3 

35 

Churchill.  Mon. 
5B.75N  94.07W 

CN 

790816 

13.52* 

80  nmi 

73.  1 
76.6 

84.  2 
87.6 

89.5 

96.0 

95.5 

100.7 

110.7 

115.6 

94.4 

98.7 

78.6 

83.9 

70.0 

73.9 

S- ' 

W.  Palm  Beach, 
26.  68N  80.  10W 

FL 

790322 

-13.63* 

63  nmi 

76.6 

81.  0 

88.4 

92.2 

94. 1 
100.6 

100. 1 
105.4 

116. 1 
120.6 

98.4 

103.3 

83.0 

88.  2 

73.3 

78.  3 

s 

W.  Palm  Beach. 
28.68N  80. 10W 

FL 

790408 

-18.  10* 

33  nmi 

74.0 

77.  1 

86.3 

89.  2 

92.0 

98.6 

98.9 

104.2 

118.3 

123.2 

96.  1 
102.9 

60.  1 
84.7 

70.6 

74.1 

>  ’  n 

Boothville,  LA 
29.  33N  89.40W 

790618 

-32.08* 

92  nmi 

81.9 

83.9 

96.3 

95.8 

104.0 

105.7 

111.6 
111.  1 

129.4 

130.6 

109.3 

106.9 

90.2 

91.3 

77.6 

80.0 

Ha-#: 

Wallops.  VA 
37.85N  75.48W 

790502 

-13.58* 

79  nmi 

76,7 

76,7 

89. 1 
89.7 

95.4 

100.0 

101.7 

105.9 

117.0 

124.9 

100.0 

104.5 

83.4 

85.0 

72.8 

73.3 

Calculated 

37.85N  75.48W 

75.0 

66.  9 

96.2 

101.5 

117.0 

99.8 

82.8 

71.9 

Measured 

Calculated 

Sterling,  VA 
38.98N  77.47W 

790518 

-13.63* 

99  nmi 

74,  1 
74.9 

86.  3 
86.7 

92.4 

96.0 

98.6 

101.4 

115.3 

116.7 

98.0 

99.7 

80.0 

82.6 

69.9 

71.7 

Wallops.  VA 

37.  85N  75.48W 

780916 

31.  91* 

23  nmi 

74.  9 

78.  5 

87.  2 

90.  9 

93.3 

100.8 

100.0 
106.  2 

115.6 

121.0 

98.4 

104.0 

61.  1 

87.  1 

70.9 

75.1 

Wallops.  VA 

37.  85N  75.48W 

790920 

-32. 14 

26  nmi 

77.2 

82.9 

89.6 

94.7 

96.2 

103.5 

102.9 

108.7 

120.  1 
126.8 

101.  1 
107.0 

84.3 

89.6 

74.  1 
80.5 

Salem.  IL 

38.  85N  88.  97W 

790516 

-31.91* 

97  nmi 

75.  2 
78.3 

87.6 

88.9 

93.  1 
96.6 

99.8 

101.2 

118.0 

116.4 

98.4 

99.6 

81.4 

84.7 

73.1 

76.0 

SSSftvr 

Chatham.  MA 
41.67N  69.97W 

790503 

-22.74* 

87  nmi 

71.3 

73.6 

81.5 

85.3 

86.2 
94.  3 

91.7 

99.5 

107.5 

115.0 

90.8 
98.  1 

76.6 

81.2 

67.8 

70.8 

Measured 

Calculated 

Chatham.  MA 
41.67N  69.98W 

790917 

9.05* 

92  ami 

78.8 

81.8 

90.  7 
93.6 

96.6 

102.6 

103.5 

107.6 

119.0 

12L5 

104.3 

105.0 

85.7 

89.8 

74.4 

78.8 

Measured 

Calculated 

Cape  Canaveral, 
28. 47N  80. 55W 

FL 

790130 

-4.52* 

47  nmi 

79.3 

74.6 

89.8 

84.6 

94.8 

92.3 

100.0 

96.7 

112.8 

110.4 

98.0 
95. 1 

85.2 

81.  1 

77.7 

72.  1 

Measured 

Calculated 

Porto  Allegro,  BZ 

30.  OS  51.  18W 

790320 

4.52* 

49  nmi 

77.7 

81.5 

91.7 

93.8 

98.  1 
103.3 

104.9 

108.5 

121.3 

125.7 

102.9 
106.  2 

84.8 

89.7 

75.0 

78.  3 

Measured 

Calculated 

Porto  Allegro,  BZ 

30.  OS  51.  18W 

790505 

27.  33* 

71  nmi 

75.0 

78.5 

89.  1 

90.  8 

96.6 

100.3 

*103.3 

105.6 

120.5 

120.6 

102.0 

103.7 

83.3 

86.7 

70.2 

75.3 

Measured 

Calculated 

Ship  Charlie 

52.  7N  35.  5W 

790507 

-8. 99* 

48  nmi 

70.3 

73.2 

81.  1 
83.3 

85.9 

90.8 

91.7 

95.3 

105.  1 
109.3 

90.9 

94.0 

76.0 

79.5 

67.0 

70.8 

Measured 

Calculated 

Ship  Charlie 

52.  7N  35.  5W 

790628 

9.  05* 

93  nmi 

72.8 

72.7 

84.4 

82.  8 

89.6 

90.3 

95.6 

94.8 

109.3 

109.6 

94.8 

93.8 

78.5 

78.7 

69.6 

70.  3 

Measured  Mean 
STD  DEV 

75.3 

3.3 

87.0 

4.3 

93.0 

4.9 

99.3 

5.4 

114.6 

7.6 

97.9 

5.2 

81.7 

4.0 

71.8 

3.3 

Calculated  Mean 
STD  DEV 

77.1 

4. 1 

88.5 

4.  8 

97. 1 
5.6 

102.  1 
6.0 

117.  3 

8.  1 

100.4 

5.6 

84.4 

4.6 

74.3 

3.9 

(Meas-Calc)  „ 
- Reas - % 

-2.4 

-1.6 

-4.4 

-2.8 

-2.3 

-2.5 

-3.3 

-3.4 

STD  DEV 

3.7 

3.4 

3.5 

3.3 

3.  1 

3.  1 

3.6 

3.9 

Correlation; 

0.697 

0.779 

0.  807 

0.834 

0.  902 

0.835 

0.758 

0.653 

i  .Tt.  -M****—^ 


Table  7.  Comparison  of  Measured  With  Calculated  Radiances— Arctic 
(62-90N)— Clear 


Location 

Date 

$ 

Distance 

FI 

355 

F2 

397 

F3 

420 

F4 

441 

F6 

535 

F6 

408 

F7 

374 

F8 

353 

Measured 

Calculated 

Thule  AB.  GL 
76.52N  68.03W 

780205 

-9.0 

18  nmi 

64.6 

68.3 

69.6 

73.5 

72.3 

76.2 

74.4 

77.9 

75.4 

77.7 

73.6 

76.6 

68.2 

71.8 

62.7 

66.6 

Measured 

Calculated 

Thule  AB.  GL 

76. 5ZN  68. 83W 

790205 

-4.5 

36  nw»i 

63.8 

68.3 

68.0 

73.5 

70.4 

76.2 

72.0 

77.9 

71.6 

77.7 

70.8 

76.6 

66.8 

71.8 

62.2 

68.8 

Measured 

Calculated 

Thule  AFB.  GL 
78.52N  88.  83W 

790205 

-13.6 

41  nml 

64.9 

68.3 

59.3 

73.5 

71.4 

76.2 

72.6 

77.9 

73.5 

77.7 

72.  1 
76.6 

67.8 

71.8 

61.0 

66.8 

Measured 

Calculated 

Khatanga.  USSR 
71.88N  102. 47E 

790418 

18.0 

75  nmi 

68.6 

70.  9 

76.9 

78.7 

80.6 

83.9 

84.8 

86.6 

88.0 

88.5 

83.5 

85.0 

73.0 

76.2 

65.3 

68.9 

Measured 

Calculated 

Khatanga.  USSR 

71.  98N  102. 47E 

790419 

13.5 

67  nmi 

69.7 

69.0 

75.7 

78.7 

78.9 

81.8 

82.0 

84.6 

84.3 

86.8 

81.1 

83.2 

72.9 

74.2 

65.9 

67.  1 

Measured 

Calculated 

Alert.  NT.  CN 

82.  5N  62.  33W 

790314 

-27.4 

23  nmi 

70  6 
70.4 

76.5 

77.2 

78.6 

81.4 

80.2 

83.2 

76.7 

81.8 

78.2 

81.6 

75.2 

75.8 

68.8 

66.4 

Measured 

Calculated 

Alert,  NT.  CN 

790321 

-22.7 

54  nmi 

66.2 

66.9 

70.9 

72.6 

73.2 

75.9 

74.6 

77.4 

72.0 

76.8 

73.6 

76.0 

69.7 

71.0 

64.3 

6S.5 

Measured 

Calculated 

Alert.  NT.  CN 

790326 

-27.3 

64  nml 

67.2 

68.  7 

74.4 

75.7 

77.7 

80.0 

80.6 

82.3 

79.7 

81.5 

78.5 

80.8 

71.5 

73.6 

65.4 

67.0 

Measured 

Calculated 

Alert.  NT.  CN 

790326 

-32.0 

98  nmi 

67.4 

68.7 

74.7 

75.7 

78.0 

80.0 

80.7 

82.3 

80.6 

81.5 

79.6 

80.8 

71.3 

73.6 

65.8 

67.0 

Measured  Mean; 
STD  DEV 

66.9 

2.2 

72.9 

3.4 

75.7 

3.8 

78.0 

4.6 

78.0 

5.8 

76.9 

4.5 

70.7 

2.8 

64.7 

2.2 

Calculated  Mean: 
STD  DEV 

68.8 

1.  2 

75.2 

2. 1 

79.  1 

3.0 

81.  1 

3.4 

81.  1 

4.2 

79.6 

3.3 

73.3 

1.9 

67.  1 

1.0 

(Meas-Calc)  m 

-2.8 

Hili  * 

Correlation: 

0.  688 

0.919 

0.961 

0.969 

0.961 

0.  964 

0.  881' 

0.524 

i 

l 
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Figure  4.  Measured  and  Calculated  Clear  Column  Radiances  for 
DMSP  SSH  H00  Channels— Tropical 
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Figure  5.  Measured  and  Calculated  Clear  Column  Radiances  for 
DMSP  SSH  H20  Channels— Mid  -  latitude 
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Figure  6.  Measured  and  Calculated  Clear  Column  Radiances  for 
DMSP  SSH  Ho0  Channels— Arctic 
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6.2  Cloud  Contaminated  Water  Vapor  Radiance  Companion 

The  colocated  cloud  contaminated  cases  selected  for  comparison  are  listed  in 
Table  8  (Tropical),  Table  9  (Mid-latitude— Low  Clouds),  Table  10  (Mid-latitude— 

High  overcast)  and  Table  11  (Arctic).  The  29  cloud  contaminated  cases— Tropical 
(5),  Mid-latitude-Low  Clouds  (11),  Mid-latitude-High  overcast  (7),  and  Arctic  (6) 
were  selected  based  again  on  the  surface  reports  and  the  3D  NEPH  data.  In  a  few 
cases  cloud  conditions  from  3D  NEPH  data  were  not  available  but  are  listed.  In 
these  tables,  the  row  designated  calculated-clear  means  that  the  transmittance  calcu¬ 
lations  were  carried  out  as  if  there  were  no  clouds  present.  In  the  case  where  the 
transmittance  calculations  were  computed  and  an  effective  cloud  layer  was  deter¬ 
mined  in  the  row  designated  calculated- cloud  at  some  height,  the  following  assump¬ 
tions  were  made:  (1)  the  cloud  has  an  emissivity  of  1.  00  and  is  constant  over  the 
width  of  the  filter  function,  (2)  the  cloud  amount  is  100  percent  and  fills  the  field  of 

view  of  the  instrument,  (3)  T  ,  the  surface  temperature  in  Eq.  (1)  is  replaced  by 

s 

T  ,  the  cloud  top  temperature,  and  (4)  the  radiance  percentage  deviation  is  mini- 
c  _  i 

mized  by  iterations  for  DMSP  channel  F5-535  cm  only,  and  (5)  the  effective  cloud 

level  found  for  F5-535  cm-1  is  valid  for  all  the  DMSP  water  vapor  channels. 

Figures  7,  8,  9,  and  10  show  the  results  of  the  comparison  in  graphical  form 
for  the  calculated  clear-cloud  contaminated  for  Tropical,  Mid-latitude-Low  Clouds, 
Mid-latitude— High  overcast  and  Arctic,  respectively.  As  should  be  expected,  in 
the  majority  of  cases,  the  calculated  exceeds  the  satellite  measured  water  vapor 
radiances.  The  higher  the  effective  cloud  layer,  the  greater  the  deviation.  This 
should  be  expected  since  the  upwelling  radiation  is  coming  from  a  higher  level  in 
the  atmosphere.  However,  it  is  surprising  that  even  with  the  cloud  contamination 
there  are  a  few  comparisons  showing  a  slight  positive  percentage  deviation.  For 
example  in  Table  8  and  Figure  9,  the  comparison  radiance  data  for  Bangkok  would 
indicate  a  clear  cdlumn  condition  but  the  3D  NEPH  data  showed  a  broken  cirrus 
stratus  condition  1.5  to  4.  0  km  thick.  No  effective  cloud  layer  calculations  were 
performed  for  Bangkok  because  the  radiance  data  for  F5-535  cm’ 1  was  already 
maximized  with  minimum  percentage  deviation.  Also  shown  in  Table  8,  the  greatest 
discrepancy  found  in  this  study  was  Ponope,  Caroline  Islands.  The  range  of  dis¬ 
crepancies  range  from  -61  percent  for  F8-353  cm  the  most  opaque  channel,  to 
-182  percent  for  F5-535  cm  *,  the  most  transparent  channel.  Although  there  are 
only  five  cases  in  the  Tropics  used  for  comparison,  the  cloud  retrievals  are  dis¬ 
appointing  when  the  effective  cloud  layer  is  compared  directly  to  the  cloud  condi¬ 
tions  reported  in  the  3D  NEPH  data.  For  example,  in  Table  8  see  Bangkok  on  790403 
and  either  Truk  or  Ponape.  The  effective  cloud  layer  calculations  do  not  even  come 
close  to  the  actual  cloud  conditions  being  reported  by  the  3D  NEPH  data.  However, 

I  the  relationship  between  the  measured  and  calculated  radiances  for  an  effective 

1 

cloud  layer  in  the  Tropica  ia  greatly  improved  a8  §hown  graphically  in  Figure  11* 

I 

i 


24 


Comparison  of  Measured  With  Calculated  Radiances— Mid-latitude  (26-62N)  Calculated  Clear  and  Effective 
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Table  11.  Comparison  of  Measured  With  Calculated  Radiances— Arctic  (62-90N) 
Calculated  Clear  and  Clouds 


FI 

F2 

F3 

F4 

F  5 

F6 

F7 

F8 

Location 

Date 

0 

Distance 

3S5 

397 

420 

441 

535 

408 

374 

353 

Measured 

Thule  AB.  GL 

790204 

0» 

5  nmi 

61.6 

63.4 

63.  2 

62.2 

60.6 

63.  2 

62.8 

60.6 

Calc-Clear 

76.  52N  68.  83W 

(13536) 

66.  e 

72.6 

75.9 

78.2 

79.3 

77.0 

70.  6 

64.  6 

Calc-Cloud  at  4.  5  km 

60.7 

63.3 

64.0 

64.  2 

61.  1 

63.7 

62.  1 

60.5 

Measured 

Thule  AB.  GL 

790294 

-4.  5® 

27  nmi 

61.  9 

64.2 

64.  1 

63.5 

62.  1 

64.3 

63.6 

61.2 

Calc-Clear 

(13536) 

66.  8 

72.6 

75.9 

78.  2 

79.  3 

77.0 

70.6 

64.  9 

Calc-Cloud  at  4.  0  km 

61.6 

64.5 

65.5 

65.9 

63.  2 

65.3 

63.3 

61.3 

Measured 

Thule  AB.  GL 

790204 

+4.5° 

37  nmi 

63.0 

66.  2 

66.9 

67.6 

67.  1 

67.3 

64.7 

62.6 

Calc-Clear 

(13536) 

66.8 

72.6 

75.9 

78.2 

79.  3 

77.0 

70.6 

64.9 

Calc-Cloud  at  3.  0  km 

63.3 

66.8 

68.4 

69.  2 

67.  3 

68.4 

65.5 

62.7 

Measured 

Thule  AB.  GL 

790204 

0° 

16  nmi 

61.  6 

63.6 

63.9 

63.0 

61.4 

63.6 

$3.5 

60.  7 

Calc-Clear 

(14537) 

66.4 

72.5 

76.  2 

78.  2 

79.  2 

77.0 

70.  5 

65.0 

Calc-Cloud  at  4.  25  km 

61.  0 

63.6 

64.4 

64.  7 

61.6 

64.2 

62.5 

60.  8 

Measured 

Thule  AB.  GL 

792004 

+4.5* 

16  nmi 

63.2 

66.6 

67.5 

67.6 

67.  1 

67.  2 

65.4 

62.5 

Calc-Clear 

(14537) 

66.4 

72.5 

76.  2 

78.  2 

79.  2 

77.  0 

70.5 

65.  0 

Calc-Cloud  at  2.5  km 

63.5 

67.  2 

68.  9 

69.6 

67.6 

68.9 

69.  8 

62.  9 

Measured 

Thule  AB.  GL 

790204 

-4.  5° 

50  nmi 

62.  0 

63.6 

64.  3 

63.7 

61.  9 

64.3 

63.  I 

60.  7 

Calc-Clear 

(14537) 

66.4 

72.5 

76.  2 

78.  2 

79.  2 

77.0 

70.5 

65.  0 

Calc-Cloud  at  4.  25  km 

61.  0 

63.6 

64.4 

64.  7 

61.6 

64.  2 

62.  5 

60.  8 

Mean- Measured 

62.25 

64.  6 

65.0 

64.6 

63.4 

65.0 

63.  85 

61.4 

-Calc -Clear 

66.6 

72.  55 

76.05 

78.  2 

79.  25 

77.0 

70.55 

64.95 

-Calc-Clouds 

61.85 

64.  8 

65.9 

66.4 

63.7 

65.8 

63.6 

61.5 

(Meas-Calc)  « 

Meas  Clear  ' 

-7.  0 

-12.  3 

-17.0 

-21.  1 

-25.  0 

-18.  5 

-10.  5 

-5.  8 

(Meas-Calc)  * 

Meas  Clouds  '* 

+  0.6 

-0.3 

-1.4 

-2.  8 

-0.5 

-1.  2 

+  0.  4 

-0.  2 

TROPICAL  (0-291) 

CALQIATED  CLEAR  -  OUT  dKOWNATB) 


«  F1-S5W 
A  R-W 
+  rc-«o 
°  BMW 
x  F5-S5 

■  FW08 
*  17-S* 

■  FJ-55 
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Figure  8.  Measured  and  Calculated  Clear  Column— Cloud 
Contaminated  Low  Cloud  Radiances  for  DMSP  SSH  H„0 
Channels— Mid-latitude 


Figure  9.  Measured  and  Calculated  Clear  Column— Cloud 
Contaminated  High  Overcast  Radiances  for  DMSP  SSH  H«0 
Channels— Mid-latitude 
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Figure  10.  Measured  and  Calculated  Clear  Column— Cloud 
Contaminated  Radiances  for  DMSP  SSH  H20  Channels— Arctic 


Figure  11.  Measured  and  Calculated  High  Overcast  Radiances  for 
DMSP  SSH  HgO  Channels— Tropical 
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The  comparison  of  measured  with  calculated  radiances  in  the  cloud  contaminated 
cases  for  the  Mid-latitude  belt  was  sub-divided  into  a  low  and  high  cloud  atmos¬ 
pheric  condition.  The  low  cloud  category  was  arbitrarily  chosen  based  on  the  cal¬ 
culation  of  an  effective  cloud  layer  at  or  below  4  km  or  approximately  the  600  mb 
level  in  the  atmosphere.  The  actual  cloud  conditions  both  in  amounts  and  levels  as 
obtained  from  the  3D  NEPH  data  are  quite  variable  as  can  be  seen  in  Table  9.  The 
Mid-latitude— High  overcast  category  shows  a  more  uniform  type  of  cloud  conditions, 
as  seen  in  Table  10. 

The  Mid-latitude  — Low -effective  cloud  layer  comparison  is  listed  in  Table  9  and 
shown  graphically  in  Figure  8.  The  mean  percentage  deviations  are  negative  for 
all  DMSP  channels.  The  range  is  -3.  7  percent  to  -5.  7  percent  for  the  cloud  con¬ 
taminated  cases  which  is  very  similar  to  the  range  of  -1.6  percent  to  -4.4  percent 
for  the  clear  column  comparison  shown  in  Table  6.  Also,  it  would  be  very  difficult 
to  see  any  significant  differences  between  Figure  5— Mid-latitude  Clear  and  Fig¬ 
ure  8— Mid-latitude  Calculated  Clear-Cloud  contaminated.  The  cloud  retrievals  only 
show  a  slight  improvement  as  shown  in  Table  9  and  Figure  12.  The  maximum  im¬ 
provement  is  shown  in  DMSP  channel  F5-535  cm  *,  a  reduction  from  -5.  7  to  -0.  3 
percent,  and  this  improvement  is  due  to  the  design  of  the  cloud  retrieval  calcula¬ 
tion.  There  is  very  little  improvement  in  the  more  opaque  water  vapor  channels, 
that  is,  F8-353  cm  1  and  Fl-355  cm  *.  It  appears  from  this  analysis  that  the  DMSP 
water  vapor  channels  cannot  distinguish  between  clear  and  low  cloud  contamination. 
Low  cloud  contamination  being  an  overcast  below  3  or  4  km  or  lower  than  the  600  to 
700  mb  atmospheric  level. 

On  the  other  hand,  the  Mid-latitude —High  overcast  comparison  listed  in  Table  10 
rn d  shown  graphically  in  Figure  9  shows  a  definite  negative  discrepancy  between 
the  measured  and  calculated  water  vapor  radiances.  The  range  of  negative  devia¬ 
tions  is  -13.0  percent  for  F8-353  cm  1  to  -61.6  percent  for  F5-535  cm  .  As 
should  be  expected,  there  is  a  definite  negative  discrepancy  as  shown  both  in 
Tables  10  and  Figure  9.  The  cloud  retrievals  show  an  improvement  in  the  compari¬ 
son  of  effective  cloud  layer  to  the  cloud  conditions  reported  in  the  3D  NEPH  data 
for  the  Mid -latitude— High  overcast  category.  The  cloud  retrievals  calculations  in 
the  measured  and  calculated  radiance  comparison  show  a  large  improvement  as 
shown  in  Table  10  and  Figure  13.  The  negative  deviation  for  channel  F5-535  cm’1 
is  reduced  from  -61.6  percent  to  -0.  2  percent  and  for  channel  F8-353  cm  1  from 
-13.  0  percent  to  0  percent. 
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Finally,  in  Table  11  and  Figure  10,  the  comparison  of  measured  with  calcu¬ 
lated  radiances  for  the  Arctic  latitude  belt  is  shown.  There  was  no  cloud  informa¬ 
tion  available  so  no  actual  comparisons  could  be  made  between  the  calculated  effec¬ 
tive  and  the  actual  cloud  layer.  This  is  also  the  only  time  when  the  comparison 
used  both  DMSP  satellites,  that  is,  13536  and  14537.  Figure  10  shows  quite 
graphically  the  cloud  contamination  because  of  the  large  negative  deviations  of  cal¬ 
culated  being  greater  than  the  measured  water  vapor  radiances.  The  range  of  nega 
tive  deviations  are  -5.  8  percent  for  F8-353  cm" 1  and  -25.  0  percent  for  F5-535  cm 
The  Arctic  cloud  retrievals  are  shown  in  Table  11  and  Figure  14. 
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Figure  14.  Measured  and  Calculated  Overcast  Radiances  for 
DMSP  SSH  HzO  Channels  —  A retie 


7.  CONCLUSIONS 


Comparisons  have  been  made  between  DMSP  SSH  HgO  measurements  and 
classically-calculated  water  vapor  radiances  based  on  colocated  in  space  and  time 
radiosonde  observations.  Approximately  70  sets  of  radiance  data  were  compared 
both  in  clear  and  cloud  contaminated  atmospheric  conditions.  In  addition,  the  sets 
were  divided  into  three  latitude  belts,  that  is.  Tropical  (26S-26N),  Mid-latitude 
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(26-62N)  and  Arctic  (62-90N).  Systematic  discrepancies,  by  approximately  a  4:1 
ratio,  are  shown  for  all  of  the  DMSP  SSH  HgO  channels  in  comparison  to  the  cal¬ 
culated  water  vapor  radiances.  The  calculated  radiances  generally  exceed  the 
measured  radiances  with  DMSP  SSH  F3-420  cm'*  channel  showing  the  largest  dis¬ 
crepancy.  In  the  mean,  the  radiance  comparison  indicates  a  systematic  discrepancy 
less  than  5  percent  for  the  clear  column  conditions.  Results  of  this  analysis,  that 

3 

is,  calculated  exceeds  measured  radiances,  is  in  agreement  with  McClatchey's 
results  in  his  analysis  of  the  15  )im  carbon  dioxide  channels. 

In  the  clear  column  comparisons,  the  discrepancies  appear  to  be  latitudinally 
dependent.  Smaller  discrepancies  are  found  in  the  Tropics  and  the  larger  dis¬ 
crepancies  are  found  in  the  Arctic  latitude  belt.  This  may  be  due  to  the  moisture 
and  temperature  profiles  representing  the  latitude  belts,  that  is.  Tropics— more 
moisture— higher  atmospheric  temperature  and  Arctic— less  moisture  and  colder 
atmospheric  temperature.  As  would  be  expected,  the  discrepancies  are  larger  in 
the  cloud  contaminated  cases  because  the  upwelling  radiation  is  coming  from  the 
cloud  top  level.  However,  in  the  low  cloud  contamination  (less  than  600  mb)  the 
discrepancies  are  very  similar  to  the  clear  column  comparisons.  Thus  it  appears 
that  the  DMSP  SSH  HjO  channels  cannot  discriminate  between  low  cloud  contamina¬ 
tion  and  clear  column  conditions.  A  cloud  retrieval  procedure  was  used  to  improve 
the  relationship  between  the  measured  and  calculated  cloud  contaminated  water 
vapor  radiances.  However,  when  the  calculated  effective  cloud  level  is  compared 
to  the  3D  NEPH  cloud  conditions,  the  results  are  disappointing.  In  general,  the 
effective  cloud  level  is  much  higher  than  that  given  in  the  3D  NEPH  data. 

In  view  of  these  results,  it  appears  that  additional  research  should  be  done  on 
the  Forward  Problem,  that  is,  matching  measurements  with  calculations.  The 
single  point  or  scan  spot  station  comparison  from  Polar  orbiting  satellites  has  been 
disappointing.  There  is  a  definite  need  to  delineate  the  Forward  Problem  dis¬ 
crepancy.  Hopefully  a  study  of  spatial  and  temporal  variations  of  the  satellite 
radiance  data  and  carefully  selected  "ground  truth"  sites  will  answer  some  of  these 
discrepancies.  It  appears  that  a  carefully  designed  research  program  using  a  geo¬ 
stationary  satellite  sounder  such  as  the  VISSR  Atmospheric  Sounder  could  provide 
some  of  the  answers  to  the  Forward  Problem  discrepancy. 
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Appendix  A 

DMSP  SSH  H20  Filter  Trammission  Curves  and 
Digitized  Filler  Functions 


Filter  transmission  curves  and  digitized  filter  functions  for  the  eight  (8)  DMSP 

SSH  water  vapor  channels  are  provided  in  Appendix  A.  These  curves  and  the  digi- 

9 

tized  filter  functions  were  obtained  from  a  Barnes  Engineering  Company's  Report 
and  are  valid  for  Flight  II  SSH  package  aboard  spacecraft  WC  13536  launched  in 
July  1977.  These  curves  were  used  to  calculate  the  water  vapor  transmittances  and 
the  resulting  weighting  functions. 


9.  Barnes  Engineering  Co.  (1976)  Flight  II-Supplementary  Sensor  H  (SSH) 
Radiometric  Performance.  Report  Z41S-tA-6i2. 
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Band  F2,  397  Wavenumber  Channel 
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Band  F4,  441  Wavenumber  Channel  Band  F5,  535  Wavenumber  Channel 
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Band  F5.  535  Wavenumber  Channel 
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Band  F6,  408  Wavenumber  Channel 
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Band  F7,  374  Wavenumber 


